Phase diagram of a rapidly-rotating two-component Bose gas 
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We derive analytically the phase diagram of a two-component Bose gas confined in an anhar- 
monic potential, which becomes exact and universal in the limit of weak interactions and small 
anharmonicity of the trapping potential. The transitions between the different phases, which con- 
sist of vortex states of single and multiple quantization, are all continuous because of the addition 
of the second component. 
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One of the remarkable features of cold atomic gases is 
their high degree of tunability, allowing for a precise and 
flexible control over most of the experimental parame- 
ters, and paving the way for many different potential 
applications of these systems. On the experimental side, 
remarkable progress has been made by, e.g., the realiza- 
tion of confining potentials of various functional forms 
(see, e.g., [H, 0), or by the creation of mixtures of differ- 
ent atomic species, see, e.g., 

Two-component rotating Bose-Einstein condensates 
have been investigated thoroughly for the case of har- 
monic confinement, One of the most interesting phe- 
nomena is the presence of so-called "coreless vortices" 
that occur when only one of the two components car- 
ries all the angular momentum and forms a vortex state 
around the second one, which remains at rest at the vor- 
tex core 0]. Indeed, the presence of the second, non- 
rotating component gives rise to an effective anharmonic 
potential acting on the rotating component, and there- 
fore allows the formation of multiply-quantized vortex 
states [1, [H . These multiply-quantized vortex states are 
not energetically favorable in a single-component system 
in the case of harmonic confinement. 

The rotational properties of a single-component Bose- 
Einstein condensate in the presence of an anharmonic 
potential have been addressed previously, employing the 
Gross-Pitaevskii mean-field approach, or the method of 
exact diagonalization [g - 17 1 . In the case of harmonic con- 
finement the rotational frequency of the trap Q is limited 
by the trap frequency uj because of the centrifugal force: 
as O -> u the system enters a highly-correlated regime 
while for Q > uj the system is not bounded. On 
the other hand, in the case of anharmonic confinement 
the system is bounded for any value of f2. 

In this study we consider a mixture of two Bose- 
Einstein condensates, which are confined in an anhar- 
monic potential 2l| . We investigate the rotational prop- 
erties of this system as a function of the rotational fre- 
quency of the trap and the coupling between the atoms. 



As we show, in the limit of weak interactions and small 
anharmonicity of the confining potential one can derive 
the corresponding phase diagram analytically solving a 
quadratic algebraic equation. Remarkably the phase di- 
agram is exact and universal in these limits. 

Our paper is organized as follows. We first focus on 
the case of zero and sufficiently weak interatomic inter- 
actions, where the order parameters of the two species are 
multiply-quantized vortex states. The simplicity of these 
states then allows us to investigate their stability as the 
coupling constant between the atoms increases, deriving 
the phase diagram of the system as a function of the in- 
teraction strength and of the rotational frequency of the 
trap. We finally analyse and interpret our results physi- 
cally and compare them with those of a single-component 
system, as the inclusion of a second component changes 
the corresponding phase diagram rather drastically. 

Model. — We consider a mixture of two distinguishable 
bosonic atoms, labelled A and B with equal mass M, but 
with different number of atoms Na and Ng ■ The system 
is confined in a two-dimensional anharmonic potential of 
the form 



V(p) = l -Mu 2 p 2 [l 



A(^) 2 ] 
ao 



(1) 



where p is the radial coordinate in cylindrical coordi- 
nates, uj is the trap frequency, ao = y / h/(Muj) is the 
oscillator length, and A is a positive dimensionless pa- 
rameter measuring the strength of the anharmonicity of 
the trapping potential. Along the axial direction, the 
density is assumed to be homogeneous within a width Z, 
with a total density per unit length a = (Na + Nb)/Z. 
The intra- and inter-species interaction is modelled as 
a hard-core potential with scattering lengths for elastic 
atom-atom collisions claa, «bb, and a^s, which are as- 
sumed to be repulsive. The general formalism is given 
for any value of the scattering lengths, while the fi- 
nal results are presented for equal scattering lengths, 
claa = cibb = olab = a. The dimensionless parameter 
era thus gives the "strength" of the interatomic coupling. 
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FIG. 1: (Colour online) Phase diagram, where the x axis is 
the frequency of rotation of the trap Q/cj and the y axis is the 
coupling era, for equal scattering lengths. The lines show the 
discontinuous transitions between the states (m, n) = (^a = 
$m; *s = $n). Here A = 0.005 and N A /N B = 2. 



Within the mean-field approximation the two order pa- 
rameters ^>a and obey the following coupled nonlin- 
ear, Gross-Pitaevskii-like, differential equations, which in 
the rotating frame have the form 



2M 



+ V{p)*i + {g u \%\ 2 + <i,, y \' , 'M', 



Here, xa,b = Na,b/N and = a^ja. The criti- 
cal frequencies for transitions between different states 
of (to, n) can be calculated by comparing the energies 
in the rotating frame. Figure 1 shows the correspond- 
ing phase boundaries for a fixed population imbalance 
Na/Nb = 2 and a weak anharmonicity, with A = 0.005. 
In the absence of interactions the corresponding critical 
frequencies are degenerate. However, these degeneracies 
are lifted by the interactions, as seen in Fig. 1. For a 
fixed coupling and increasing 51 the system first under- 
goes a (discontinuous) transition to the state where the 
component with smaller number of atoms carries all the 
angular momentum and the other component remains 
static, which is a so-called coreless vortex state. Then, 
with increasing Q, a vortex state forms in the larger com- 
ponent, while the smaller one becomes static. Finally, a 
vortex state forms in both components for a sufficiently 
large value of ft. It is interesting that the same result 
is obtained in the case of mixtures in a purely harmonic 
potential 

Continuous transitions. — As the interaction strength 
increases, the states of multiple quantization become un- 
stable, as the energy of the system is minimized by mix- 
ing states of different angular momentum in the order 
parameters ^a and Wg. The multiply-quantized vor- 
tex states undergo continuous, second-order phase tran- 
sitions. It turns out that the order parameters above the 
phase boundaries are of the form 



(2) 



- q <f>m+q, *B = d n <f> n + d n+q <£ n+q . (4) 



with (i — A,B;j = B.A). Here L z is the axial compo- 
nent of the angular momentum operator, ^ is the chem- 
ical potential of each component and gij — Anh 2 aaij/M. 

Discontinuous transitions. — Starting with the case of 
zero coupling, a a = 0, and small anharmonicity, Ac 1, 
the order parameters ^ a and 'J b are given by the eigen- 
states of the harmonic potential with no radial nodes 
$ m (p,(f)) oc p m e lmr t> e ~p /2q 0j where m is the quantum 
number that corresponds to the angular momentum mh 
(assumed to be positive). The single-particle energy 
spectrum E m of the above states $ m scales quadratically 
with m, E m = Uu}[l+m+\(m+l)(m+2)/2], as opposed 
to the harmonic potential, where E m oc m. As a result, 
as fl increases, the system undergoes discontinuous phase 
transitions between the states <I> m . 

For weak coupling, era <C 1, the effect of interactions 
may be treated perturbatively. The energy of the system 
in a state of (to,tx), where component A is in the state 
<&,„, and component B is in the state $ n , is 
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Sufficiently close to the phase boundary, the coefficients 
c m and d n are of order unity, while the other two coeffi- 
cients in Eq. (HJ tend to zero. For this reason we keep in 
the energy only the terms which are up to quadratic in 
c-m+q and 



£ - -^{E m c 2 m + E m+q c 2 m+q ) + -^i E nd 2 n + E n+q d 2 n+q ) 
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where Sa 



^m+q 



and Sb 



d 2 n+q . Since 



the "pure" states ^> a = $m and ^b = ®n provide an 
extremum of the energy at the phase boundaries, the 
first-order derivatives of the energy in the rotating frame, 
£rot = £ — (m + n)hfl, with respect to c m + q and d n+q 
vanish at the phase boundary. Therefore, the stability 
of the states of multiple quantization is determined by 
the eigenvalues of a two by two matrix whose elements 
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consist of the second-order derivatives of the energy with 
respect to c m+q and d n+q , i.e., M M = d 2 £ rot /dc 2 n+q , 
M 2 , 2 = d 2 £ mt /dd 2 l+q , and M lj2 = d 2 £ mt /dc m+q dd n+q . 
The resulting matrix elements along the diagonal are 



Mi 



2xA(E m+q 
-V AA 1 



E m -qTiQ)+4x 2 A [2V r 
- 4x A x B [V 7 
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AB 
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* n,n,n,n] 1 L r m,n-|-(/,m,n+g 

while the off-diagonal element is M\ 5 2 

+ 



4xAXBV A n+q «i where 



(8) 



When all the eigenvalues of this matrix are positive, the 
system is stable. However, as soon as any of the eigenval- 
ues becomes negative, an instability occurs via a second- 
order and continuous phase transition. Assuming equal 
scattering lengths, we derive the following general ex- 
pression for the phase boundaries of the "pure" states 
i&A = 3>m and ^ b = for n — to, which is given by 



1 

2V r 



E m +q — E m — qfl 



in . in.m , ni 



V, 



(9) 



The parameter q in the above expression can take any 
value, provided that to + q > 0. 

The result of the calculation described above is shown 
in the phase diagram depicted in Fig. 2, which is the main 
result of this study. The phase boundaries corresponding 
to these continuous transitions (triangular regions) re- 
main in between the straight lines of discontinuous phase 
transitions of absolute energetic stability shown in Fig. 1, 
which implies that the transitions between the phases of 
multiply quantized vortex states are no longer discontin- 
uous. As long as the interaction strength has a non-zero 
value, they are continuous. The resulting phase diagram 
obtained for the two-component case in this respect is 
different from that of the one-component case for which 
the transitions between phases are found to be both con- 
tinuous and discontinuous (l3| . 

The part of the phase boundary which represents the 
instability towards the mixed states with q = — 1 are 
vertical due to the fact that the denominator of Eq. (|9]) 
becomes zero (i.e., V m ,m,m,m = V m ,m-i,m,m-i) for this 
unstable mode. The lines corresponding to the unstable 
mode with q = 1 have negative slope, since 



(era) g=i 



TO + 1 



v, 



■[l-n + A(m + 2)]. 



(10) 



m,m : m,m 



Thus, the region of stability takes the shape of a triangle 
for the pure state with m = 1. The lines with positive 



slope cutting the other triangular regions in Fig. 2 denote 
the phase boundaries of the most unstable mode for the 
states with to > 2. Accordingly, when to = 2 and to > 
6, we find that the most unstable mode corresponds to 
q = — 2. For the intermediate phases with to = 3,4, and 
5 there is a transient regime where q = —3. 

As the rotational frequency increases, we observe that 
the region of stability of the multiply-quantized vortex 
states extends further, up to the points where the lines 
with positive slope cut the triangular regions. For suf- 
ficiently large values of to, demanding that (aa) g =i — 
(<ja) q= -2 and ignoring terms of order of l/m, we find 
that f2 fa 1 + A(to + 1). In other words, for large m, the 
lines with positive slope cut the ones with negative slope 
right at the top of the triangle. 

In the obtained phase diagram, we have determined 
triple points, where three phases coexist As 
it is seen from Fig. 2, the states (2,2), (3,3), and (4,4) 
have such triple points, where the two phase boundaries 
with positive and negative slopes intersect. For exam- 
ple, in the phase (2,2) around the point where the two 
corresponding phase boundaries cut each other, there is a 
doubly-quantized vortex state, a doubly-quantized vortex 
state with a single vortex around it, and two-singly quan- 
tized vortex states. In the fast-rotating regime, > ui, 
the effective potential due to the confinement and the 
centrifugal potential has a "Mexican hat" form, which 
leads to a hole in the density of the cloud at the trap 
center. This hole appears first at (4,4), as it is seen 
in Fig. 3, where we show schematic plots of the density 
and of the phase of the two order parameters and 
for the regions above the phases (2,2), (3,3), (4,4), 
and (5,5). The instability of the phase (3,3) is towards 
($01^3)1 which has a non-vanishing density at the cen- 
ter of the cloud. On the other hand, the instability of 
the state (4, 4) is towards ($1, $4), which does vanish at 
the origin. Therefore, beyond the phase (3, 3) there is 
always a node in the density of the cloud at the center 
of the trap. From the plots in Fig. 3 we also see that 
the density minima of the one component coincide with 
the density maxima of the other component. The total 
density of the cloud thus remains as close as possible to 
axial symmetry along these phase boundaries. 

The phase diagram that we have evaluated is exact 
and universal in the limit of weak interactions and weak 
anharmonicity of the trapping potential. In the single- 
component case the phase diagram is also universal, how- 
ever it is only partly exact [13( because there are also 
discontinuous phase transitions, where the phase bound- 
aries can only be evaluated approximately. On the other 
hand, in the present problem all the phase boundaries 
shown in Fig. 2 are continuous and thus are all exact. 
The phase diagram is universal in the sense that it is in- 
variant under changes of the degree of anharmonicity A 
of the trapping potential (provided that A <C 1), under 
proper rescaling of the two axes. 
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FIG. 2: (Colour online) Phase diagram, where the x axis is 
the frequency of rotation of the trap fl/cj and the y axis is the 
coupling aa, for equal scattering lengths. All the lines show 
continuous transitions. Here (m, m) = ('J a = 3>m;^s = 
$ m ), while ($ n ,$ m ) A ' B = = c n $ n +c m $ m ;* s = 

d„$„ +d m $ m ) denotes the states of the form of Eq. Q. Here 
A = 0.005 and N A /N B = 2. 
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FIG. 3: (Color online) Density (upper) and phase (lower 
plots) of the component A (left) and B (right plots). 
These plots correspond to the "mixed" phases (3>o, <&2) A,S , 
($o,$3) A,S , (<I>i,$4) A ' S , and ($ 2 ,$5) A ' S from left to right. 



Numerical results. — We have confirmed the phase 
boundaries shown in Fig. 2 numerically for various val- 
ues of the coupling and of the rotational frequency. Since 
the transitions are continuous, we start from the phases 
of multiple quantization found in Fig. 1 for a fixed f2 
and increase the coupling. We choose the initial state 
^ A = — ®m adding some more states with 

different angular momentum with very small amplitude 
and then propagate those in imaginary time using of a 
fourth-order split-step Fourier method [22|. Below the 
phase boundaries the amplitudes of the extra components 
decay, while above the boundary one of the small am- 
plitudes increases, indicating the formation of a mixed 
state. The result of this calculation is shown as the 
crosses in the phase diagram of Fig. 2. The numerical 
results are in good agreement with the exact analytical 
solution. However, the discrepancy between the results 
of the two methods grows rapidly with increasing A, since 
our analytical approach is perturbative. 

Summary. — In summary, we have examined the phase 
diagram of a mixture of two Bose-Einstein condensed 



gases confined in an anharmonic potential under rota- 
tion, as a function of the strength of the coupling con- 
stant and of the rotational frequency of the trap. We 
have shown that it is possible to derive the correspond- 
ing phase diagram analytically, reducing the problem to 
the evaluation of the roots of an algebraic equation of 
second degree. It is also remarkable that the presence of 
a second component makes the solution of this problem 
in a sense simpler than its one-component counterpart. 

In the case of a single component, for sufficiently 
weak interactions the cloud undergoes discontinuous 
phase transitions between phases of multiple quantiza- 
tion. Here, while there are still phases of vortex states of 
multiple quantization, the transition between them takes 
place via continuous transitions. This becomes possible 
via vortex states which enter the two components suc- 
cessively from infinity, moving continuously towards the 
trap center. 

The phase diagram we have evaluated is exact for suffi- 
ciently weak interactions and for small anharmonicity of 
the trapping potential. As long as these two assumptions 
are not violated, it is also universal. In the results of our 
study we have assumed equal scattering lengths for inter- 
and intra-species collisions. In the more general case, i.e., 
when they are not equal to each other (and/or the masses 
are unequal) , the phase boundaries are no longer straight 
lines. However, the problem is still of the same level of 
difficulty. Given the simple and systematic behaviour of 
the evaluated phase diagram, it would be interesting to 
confirm these results experimentally. 
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